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Abstract 

In this paper we study the horoball packings related to the hyperbolic 

—4 

24 cell in the extended hyperbolic space H where we allow horoballs in 
different types centered at the various vertices of the 24 cell. 

We determine, introducing the notion of the generalized polyhedral den¬ 
sity function, the locally densest horoball packing arrangement and its den¬ 
sity with respect to the above regular tiling. The maximal density is ~ 
0.71645 which is equal to the known greatest ball packing density in hy¬ 
perbolic 4-space given in IT3l . 


1 Introduction 

We consider horospheres and their bodies, the horoballs. A horoball packing B of 
H is an arrangement of non-overlapping horoballs B in El . 

‘Mathematics Subject Classification 2010: 52C17, 52C22, 52B15. 
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The definition of packing density is critical in hyperbolic space as shown by 
Boroczky 0. For standard examples also see 112X1 . The most widely accepted 
notion of packing density considers the local densities of balls with respect to their 
Dirichlet-Voronoi cells (cf. 0 and JH). In order to consider horoball packings in 
H we use an extended notion of such local density. 

Let B be a horoball in packing B, and Pei be an arbitrary point. Define 
d(P, B) to be the perpendicular distance from point P to the horosphere S = did, 
where d(P, B ) is taken to be negative when P e B. The Dirichlet-Voronoi cell 
V(B,B) of a horoball B of packing B is defined as the convex body 

V{B, B) = {Pe H n |d(P, B) < d(P, B'), MB' 6 B}. 


Both B and V are of infinite volume, so the usual notion of local density is modi¬ 
fied as follows. Let Q e dW 1 denote the ideal center of B at infinity, and take its 
boundary S to be the one-point compactification of Euclidean (n — l)-space. Let 
B'f~\r) C S be an (n — l)-ball with center C G S \{Q}. Then 0 <G dW 1 and 
B(P ] (V) determine a convex cone C n (r) = coneq(B^ 1 (r)) e El with apex 0 
consisting of all hyperbolic geodesics passing through Ly'r 1 (r) with limit point 
Q. The local density S n (B, B) of B to D is defined as 


S n (B,B) 


-— vol(B D C n (r )) 
r-> 0 o vol(T> (T C n (r )) 


This limit is independent of the choice of center C for B' ( P ] (r ). 

For periodic ball or horoball packings the local density defined above can be 
extended to the entire hyperbolic space. This local density is related to the simpli- 
cial density function that was generalized in I134H and If35ll . In this paper we will 
use the generalization of this definition of packing density. 

In /[??!/ we have refined the notion of the „ congruent” horoballs in a horoball 
packing to the horoballs of the ,,same type” because the horoballs are always 
congruent in the hyperbolic space H , in general. 

Two horoballs in a horoball packing are in the „same type”, or ,,equipacked”, 
if and only if the local densities of the horoballs to the corresponding cell (e.g. 
D-V cell; or ideal regular polytop, later on) are equal. 

If we assume that the „horoballs belong to the same type”, then by an¬ 
alytical continuation, the well known simplicial density function on H can be 
extended from n-balls of radius r to the case r = oo, too. Namely, in this case 
consider n + 1 horoballs which are mutually tangent and let B be one of them. 
The convex hull of their base points at infinity will be a totally asymptotic or ideal 
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regular simplex T™ g 
to write 


G H of finite volume. Hence, in this case it is legitimated 


d n ( OO) 


(n + 1) 


vol(B n T r ”) 

vol(T™ g ) 


Then for a horoball packing B, there is an analogue of ball packing, namely (cf. 
||4j| . Theorem 4) 

5 n (B,B)<d n ( oo), VU e B. 


Remark 1.1 The upper bound d n ( oo) (n = 2, 3) is attained for a regular horoball 
packing, that is, a packing by horoballs which are inscribed in the cells of a reg¬ 
ular honeycomb of H . For dimensions n = 2, there is only one such packing. It 
belongs to the regular tessellation {oo, 3} . Its dual {3, oo} is the regular tessel¬ 
lation by ideal triangles all of whose vertices are surrounded by infinitely many 
triangles. This packing has in-circle density d 2 (oo) — - m 0.95493. 

_ 3 

In H there is exactly one horoball packing with horoballs in same type whose 
Dirichlet-Voronoi cells give rise to a regular honeycomb described by the Schlafli 
symbol {6, 3, 3} . Its dual {3, 3, 6} consists of ideal regular simplices Tff g with 
dihedral angle | building up a 6-cycle around each edge of the tessellation. The 
density of this packing is 6f ~ 0.85328 

If horoballs of different types at the various ideal vertices are allowed 

i.e the horoballs are differently packed, then we generalized the notion of the 
simplicial density function lf34ll . In lfl2ll we proved that the optimal ball packing 
arrangement in HI 3 mentioned above is not unique. We gave several new examples 
of horoball packing arrangements based on totally asymptotic Coxeter tilings that 
yield the Boroczky-Florian upper bound If5il . 

Furthermore, in |[34]| . 051 we found that by admitting horoballs of different 
types at each vertex of a totally asymptotic simplex and generalizing the simplicial 
density function to H for [n > 2), the Boroczky-type density upper bound is no 

- 4 

longer valid for the fully asymptotic simplices for n > 3. For example, in El the 
locally optimal packing density was found to be 0.77038 ... which is higher than 
the Boroczky-type density upper bound 0.73046 .... However these ball packing 
configurations are only locally optimal and cannot be extended to the entirety of 
the hyperbolic spaces EI n . 

In [fl3l we have continued our investigations on ball packings in hyperbolic 
4-space. Using horoball packings, allowing horoballs of different types, we find 
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seven counterexamples with density « 0.71645 (which are realized by allowing 
up to three horoball types) to one of L. Fejes-Toth’s conjectures. 

Several extremal properties relate to the regular hyperbolic 24-cell and the 
corresponding Coxeter honeycomb concerning the right angled polytops and hy¬ 
perbolic 4-manifolds. 

A. Kolpakov in imi has shown that the hyperbolic 24-cell has minimal volume 

—4 

and minimal facet number among all ideal right-angled polytopes in H . 

J. G. Ratcliffe and S. T. Tschantz in ll22l have constructed complete, open, hy¬ 
perbolic 4-manifolds of smallest volume by gluing together the sides of a regular 
ideal 24-cell in hyperbolic 4-space. They also showed that the volume spectrum 
of hyperbolic 4-manifolds is the set of all positive integral multiples of Arc 2 /3. 

L. Slavich has constructed in Il24l . using the hyperbolic 24-cell, two new ex¬ 
amples of non-orientable, noncompact, hyperbolic 4-manifolds. The first has min¬ 
imal volume V m = 47r 2 /3 and two cusps. This example has the lowest number 
of cusps among known minimal volume hyperbolic 4-manifolds. The second has 
volume 2 • V m and one cusp. It has lowest volume among known one-cusped 
hyperbolic 4-manifolds. 

In This paper we study a new extremal property of the hyperbole regular 24- 
cell and the corresponding regular 4-dimensional honeycomb described by the 
Schldfli symbol {3, 4, 3, 4} relating to horoball packings. 

We determine, introducing the notion of the generalized polyhedral density 
function, the locally densest horoball packing arrangements and their densities 
with respect to the above 4-dimensional regular tiling. The maximal density is 
~ 0.71645 which is equal to the known greatest ball packing density in hyperbolic 
4-space given in m 


2 Formulas in the projective model 

We use the projective model in Lorentzian [n + l)-space E 1,n of signature (1, n), 
i.e. E 1,n is the real vector space V n+1 equipped with the bilinear form of signature 

(l,n) 

( x, y) = -x°y° + x l y 1 4- b x n y n (2.1) 

where the non-zero real vectors x = (x°, x 1 ,..., x n ) £ V n+1 and y = ( y°, y 1 ,..., 
y n ) e V n+1 represent points in projective space V n (W). 1" is represented as the 
interior of the absolute quadratic form 


Q = {[x] G P n |(x, x) = 0} = dW 1 


(2.2) 
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in real projective space V n (V n+1 , V n+ \). All proper interior points x G H n are 
characterized by ( x, x) < 0. 

The boundary points dW n in V n represent the absolute points at infinity of HI"'. 
Points y satisfying ( y, y) > 0 lie outside dW 1 and are called the outer points of 
H”. Take P([x]) G V n , point [y] G V" is said to be conjugate to [x] relative to 
Q when ( x, y) = 0. The set of all points conjugate to P([x]) form a projective 
(polar) hyperplane 


pol(P):={[ y]er|(x, y) = o}. (2.3) 

Hence the bilinear form Q in (12.11) induces a bijection or linear polarity V n+1 —> 
V n+ \ between the points of V" and its hyperplanes. Point X[x\ and hyperplane 
a[a\ are incident if the value of linear form a evaluated on vector x is zero, i.e. 
xa = 0 where x G V n+1 \ {0}, and a G V n+ \ \ {0}. Similarly, lines in V n are 
characterized by 2-subspaces of V” +1 or (n — l)-spaces of V n+ i [?]. 

Let P C HP denote a polyhedron bounded by a finite set of hyperplanes H l 
with unit normal vectors b l G V n +\ directed towards the interior of P : 

H l ■= {x G M d \( x, V) = 0} with (6\ b l ) = 1. (2.4) 


In this paper P is assumed to be an acute-angled polyhedron with proper or ideal 
vertices. The Grammian matrix G(P ) := (( b\ b 3 ))i 1 j i,j G {0,1, 2 ... n} is an 
indecomposable symmetric matrix of signature (l,n) with entries ( b 1 . b‘) = 1 
and (b\ b‘) < 0 for i =4 j where 


(b\ V) 


' 0 

— cos aP 

- 1 

! — cosh r j 


if PL i JL H j , 

if P[\ intersect along an edge of P at angle o lJ . 

if H\ Hi are parallel in the hyperbolic sense, 

if H\ II 1 admit a common perpendicular of length IP. 


This is visualized using the weighted graph or scheme of the polytope The 

graph nodes correspond to the hyperplanes H l and are connected if W and II > not 
perpendicular (/' ^ j). If they are connected we write the positive weight k where 
oiij = n/k on the edge, and unlabeled edges denote an angle of 7 t/3. 

In this paper we set the sectional curvature of H n , K = —k 2 , to b e k = 1. The 
distance d of two proper points [x] and [y] is calculated by the formula 


-(x, y) 

\/( x > x )(y> y) 


cosh d 


(2.5) 
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The perpendicular foot Y[ y] of point X[x] dropped onto plane \u\ is given by 


where [u] is the pole of the plane [u]. 

A horosphere in HP (n > 2) is a hyperbolic n-sphere with infinite radius 
centered at an ideal point on dW L . Equivalently, a horosphere is an (n - 1)- 
surface orthogonal to the set of parallel straight lines passing through a point of the 
absolute quadratic surface. A horoball is a horosphere together with its interior. 

We consider the usual Beltrami-Cayley-Klein ball model of HP centered at 
0(1, 0, 0, ..., 0) with a given vector basis a* (i = 0,1, 2,..., n ) and set an arbi¬ 
trary point at infinity to lie at T 0 = (1, 0,..., 0,1). The equation of a horosphere 
with center T 0 = (1, 0,..., 1) passing through point S = (1, 0,..., s) is derived 
from the equation of the the absolute sphere — x°x°+x 1 x 1 +x 2 x 2 + - ■ ■+x n x n = 0, 
and the plane x° — x n = 0 tangent to the absolute sphere at T 0 . The general equa¬ 
tion of the horosphere is in projective coordinates (s ±1): 



(2.7) 


and in cartesian coordinates setting h % = ftr it becomes 



( 2 . 8 ) 


1 - s (1 - s) 2 


In n-dimcnsional hyperbolic space any two horoballs are congruent in the clas¬ 
sical sense. However, it is often useful to distinguish between certain horoballs 
of a packing. We use the notion of horoball type with respect to the packing as 
introduced in [[Ml . 

Two horoballs of a horoball packing are said to be of the same type or equi- 
packed if and only if their local packing densities with respect to a given cell (in 
our case hyperbolic 24 cells) are equal. If this is not the case, then we say the two 
horoballs are of different type. 

In order to compute volumes of horoball pieces, we use Janos Bolyai’s classi¬ 
cal formulas from the mid 19-th century: 

1. The hyperbolic length L(x) of a horospheric arc that belongs to a chord 
segment of length x is 



(2.9) 
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2. The intrinsic geometry of a horosphere is Euclidean, so the (n—1) -dimensional 
volume A of a polyhedron A on the surface of the horosphere can be calcu¬ 
lated as in E n_1 . The volume of the horoball piece Hi A) determined by A 
and the aggregate of axes drawn from A to the center of the horoball is 

Vol(U(A )) = —A (2.10) 

n — 1 

3 On hyperbolic 24 cell 

An //-dimensional honeycomb V, also referred to as a solid tessellation or tiling, 
is an infinite collection of congruent polyhedra (polytopes) that fit together face- 
to-face to fill the entire geometric space (at present HP (d A 2)) exactly once. 
We take the cells to be congruent regular polyhedra. A honeycomb with cells con¬ 
gruent to a given regular polyhedron P exists if and only if the dihedral angle of P 
is a submultiple of 2n (in the hyperbolic plane zero angles are also permissible). 

A complete classification of honeycombs with bounded cells was first given by 
Schlegel in 1883. The classification was completed by including the polyhedra 
with unbounded cells, namely the fully asymptotic ones by Coxeter in 1954 [!6l. 
Such honeycombs (Coxeter tilings) exist only for d < 5 in hyperbolic //-space W l . 

An alternative approach to describing honeycombs involves analysis of their 
symmetry groups. If V is a Coxeter honeycomb, then any rigid motion moving 
one cell into another maps the entire honeycomb onto itself. The symmetry group 
of a honeycomb is denoted by SymV. The characteristic simplex T of any cell 
P E V is a fundamental domain of the symmetry group SymV generated by 
reflections in its facets which are (d — 1)-dimensional hyperfaces. 

The scheme of a regular polytope P is a weighted graph (diagram) characteriz¬ 
ing P C W l up to congruence. The nodes of the scheme, numbered by 0.1,, d, 
correspond to the bounding hyperplanes of T. Two nodes are joined by an edge if 
the corresponding hyperplanes are non-orthogonal. Let the set of weights (n 1; n 2 , 
n 3 ,..., ?T.d_i) be the Schlafli symbol of P, and n ( i be the weight describing the 
dihedral angle of P, such that the dihedral angle is equal to —. In this case T is 
the Coxeter simplex with the scheme: 

n, n 2 n d _, n d 

O-O-O . o-o-o . 

0 12 d-2 d-1 d 

Figure 1: Coxeter-Schlafli simplex scheme 
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The Schlafli symbol of the honeycomb V is the ordered set (rii, n 2 , n 3 ,..., 
rid-i,n d ) above. A (d + 1) x (d + 1) symmetric matrix (Jf 3 ) is constructed for 
each scheme in the following manner: b n = 1 and if i ^ j e {0,1,2,..., d} then 
b 13 = — cos — . For all angles between the facets i,j of T holds then n*, = rik-i,k- 
Reversing the numbering of the nodes of scheme V while keeping the weights, 
leads to the scheme of the dual honeycomb V* whose symmetry group coincides 
with SymV. 

If SymV denotes the symmetry group of a honeycomb then one tile P of 
the Coxeter tiling V nin2 ... nd can be derived by the above symmetry group and its 
characteristic simplex V: 

Pmri2—n d = | 'y(J 7 n 1 n 2 —n d ) ^ • 

7 e Sym'Pn 1 n 2 ...n d _ 1 


Every //-dimensional totally asymptotic regular polytope P has a hyperbolic 
ideal presentation obtained by normalising the coordinates of its vertices so that 
they lie on the unit sphere S n_1 and by interpreting 1 as the ideal boundary of 
H in Beltrami-Cayley-Klein’s ball model. Therefore the ideal regular hyperbolic 
24-cell P 24 can be derived from the Euclidean 24-cell as the convex hull of the 
points 


M1, 72 ,0, 72 ,0); 

A5(1, 7! ,0,0, 71 ); 

A 6 (i,- — 

A(1,0, 72 , 72 ,0); 

^(i, 0 ,—, 0 ,—); 


Al4(1, “7f’7! ,0,0); 

^(i,——,0); 

A 18 (l,—,0,0,—); 


(3.1) 
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Alo(1,0,- 7I’ 0, 7I ); 


A 24 (l,0,0,—,- — y, 


where the points (vertices) are described in a projective coordinate system given 
in Section 1. 

The 24-cell is the unique regular four-dimensional polytope having cubical 
vertex figure because the vertex figure of the other five regular four-dimensional 
polytopes are other Platonic solids, and therefore their dihedral angles are not 
sub-multiples of n thus only the regular 24-cell may be used as a building block 
in order to construct cusped hyperbolic 4-manifolds. 


3.1 The structure of the hyperbolic 24-cell P 24 

P 2 4 is a tile of the 4-dimensional regular honeycomb P 24 with Schlafli symbol 
{3,4, 3,4}. It has 24 octahedral facets, 96 triangular faces, 96 edges and 24 cu¬ 
bical vertex figures. A hyperbolic 24 cell contain 24-48 = 1152 characteris¬ 
tic simplex P 24 an d the volume of such a Coxeter simplex with Schlafli symbol 
{3,4, 3,4} is V ol(fFi 4 ) = (see [HI) therefore the volume of the hyperbolic 
24-cell is Vol(V 2 4 ) = fvr 2 . 

The vertices of P 24 are denoted by (i G {1, 2,..., 24}) and they coordinates 
are given in (3.1). 

We introduce the notion of the ^-neighbouring points (k G {1, 2, 3,4}) related 
to the vertices of P 24 : 

Definition 3.1 1. The 1-neighbouring vertices of Ai {i G {1,2,..., 24}) among 

the vertices of P 24 are the vertices Aj where AiAj is an edge o/P 24 . 

2. The 2-neighbouring vertices of Ai (1 G {1, 2,..., 24}) among the vertices 
of P 24 are Aj where AiAj is a diagonal of an octahedral facet of P 24 . 

3. The 4-neighbouring vertex of Ai (i G {1,2,..., 24}) among the vertices of 
P 24 is it opposite vertex Aj regarding P 24 . 

4. The 3-neighbouring vertices of Ai (i G {1, 2,..., 24}) among the vertices 
of P -24 are the vertices Aj that are not k-neighbouring (k = 1, 2,4) vertices 
of A^ 
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A, 



A 3 A 5 A 7 A 9 A I6 A IS A 30 A 22 


1-neighbouring vertices of A t 


A 2 A n A 12 A I4 A 23 A 24 


2-neighbouring vertices of A 7 


A 4 A 6 A s A 10 A 1s A I7 A 19 A 21 

A, 


l 13 


3- neighbouring vertices of A t 

4- neighbouring vertex of A 7 


Figure 2: The ’’neighborhood structure” of Po 4 


The Fig. 2 shows the ^-neighbouring vertices (k = 1, 2, 3,4) of A 2 . 

Definition 3.2 Two horoballs B, and Bj (or horn spheres Bf and B? ( i,j e 
{1,2,..., 24}, i 7 j ) among the horoballs centered at the vertices Po4 are k- 
neighbouring (k 6 {1,2, 3,4}) if their centres Ai and Aj are k-neighbouring 
vertices regarding Po 4 . 

We choose a characteristic simplex (orthoscheme) of P 2 4 with vertices T 0 = 
Ai (■*■> 73 ’ 73 ’ 0’ oj, T\, T 2 , T 3 and P 4 = O where P 4 ( 1, 0, 0, 0, 0) is the centre of 

P 2 4 (coincides with the center of the model), T 3 ^ 1 , A^, A—, A^, A^ j is the cen¬ 
tre of the facet-polyhedron AiA 3 A 5 A 7 A 9 An (octahedron), the centre of its regu¬ 
lar face-polygon A X A 3 A 7 (regular triangle) is denoted by T 2 (l, A^, A^, A-, 0 j 

and Ti ^1, 4j, A^, A^, (ij is the centre of the edge A t A :} of this face. Moreover, 

we denote by T^l, A^, A^, 0^ the center of the edge A 3 A 7 . This point is 

coincide with the orthogonal projection of A 2 onto its adjacent octahedral facet 

A^A^A-jA^AuA^ (see Fig. 3). 


4 Horoball packings and polyhedral density func¬ 
tion 

Similarly to the above section let P 24 be a tile of the 4-dimensional regular honey¬ 
comb P 24 with Schlafli symbol {3,4, 3,4}. We study the horoball packings B with 
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to 


o = t 4 


o 


Figure 3: A part of P 24 


horoballs centred at the infinite vertices of V 24 . The horospheres and horoball cen¬ 
tred at the vertex Ai are denoted by Bf and Bi. The density 8 (B) of the horoball 
packing B relating to the above Coxeter tiling can be defined as the extension of 
the local density related to the poly top P 2 4 - It is well known that for periodic ball 
or horoball packings the local density can be extended to the entire hyperbolic 
space. 

Definition 4.1 We consider the polytop P 24 with vertices A, (?’ = 1,.... 24) in 4- 

- 4 

dimensional hyperbolic space H . Centres of horoballs lie at vertices ofP 2 4 . We 
allow horoballs ( B^ , z = 0,1,2,..., 24) of different types at the various vertices 
and require to form a packing, moreover we assume that 



where the hyperplanes O, (j = 1,..., 18) do not contain the vertex A r . The 
generalized polyhedral density function for the above polytop and horoballs is 
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defined as 

„ R i 

Vol(P 2 4 ) 

The aim of this section is to determine the optimal packing arrangements B opt and 

- 4 

their densities for the regular honeycomb V 2 a in H . We vary the types of the 
horoballs so that they satisfy our constraints of non-overlap. The packing density 
is obtained by the above definition. 

We will use the consequences of the following Lemma (see ll35l l: 

Lemma 4.2 Let If and B 2 denote two horoballs with ideal centers Cf and C 2 , 
respectively, in the n-dimensional hyperbolic space (n > 2). Take T\ and r 2 
to be two congruent n-dimensional convex piramid-like regions, with vertices Cj 
and C 2 . Assume that these horoballs B A ( x ) and Bo ( x ) are tangent at point I ( x ) G 
C' 16'2 and CfC 2 is a common edge o/ti and r 2 . We define the point of contact I (0) 
(the so-called,,midpoint” )such that the following equality holds for the volumes 
ofhoroball sectors: 

L(0) := 2vol(B\(0 ) (T rf) = 2uo/(S 2 (0) flr 2 ). 

If x denotes the hyperbolic distance between 1(0) and l(x), then the function 

V(x) := vol(Bi(x) D Ti) + vol(B 2 (x) D r 2 ) = - (e^ n-1 ^ + e _< - n_1 ^) 
strictly increases as x — * ±cxd. 

We consider the following four basic horoball configurations £>j, (1 = 0,1, 2, 3,4): 

1. All 24 horoballs are of the same type and the adjacent horoballs touch each 
other at the ,,midpoints” of each edge. This horoball arrangement is denoted 

by B 0 . 

2. We allow horoballs of different types and the opposite horoballs e.g. Bi and 
£>13 touch their common 2-neighbouring horoballs Bi (i = 2,11,12,14, 23, 
24) (see Fig. 2) at the centres of the corresponding octahedral facets e.g. the 
horoball /J, touches the horoball If , at the facet center T 3 and f ? 13 tangent 
B u at the centre of octahedral facet A a A ( ,A h A H] A } , /I , :3 (see Fig. 3). The 
other ’’smaller” horoballs are in the same type regarding P 24 and touch their 
1 -neighbouring ’’larger” horoballs e.g. the ’’larger” horoballs B\ and If \ 
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touch the ’’smaller” horoballs P 3 , P 5 , B~. P 9 . At this horoball arrangement 
let the point A 3 A 3 n B{ be denoted by C — I\ (see Fig. 4.a) (P* is the 
corresponding horosphere of horoball P,.) 

This horoball arrangement is denoted by B\. 

3. We set out from the B\ ball configuration and we expand the horoballs B\ 
and £>i 3 until they comes into contact with their adjacent facets regarding 
P 2 4 while keeping their 1 and 2-neighbouring horoballs tangent to them. At 
this configuration which is denoted by P 2 the horoballs are included on 3 
classes related to P 24 . The horoballs Pi and Pi 3 are in the same type and 
they touch their corresponding 1 -neighbouring horoballs that form the sec¬ 
ond class. The remaining 8 horoballs are also in same type and are included 
on the 3. type. 

For example the horoball Pi touches its neighbouring facet at the point 
T (see Fig. 3, and Fig. 4.b) and touches its 1-neighbouring horoballs e.g. 
P 3 , P 5 , P 7 , P 9 and its 2-neighbouring horoballs e.g. P 14 . At this horoball 
arrangement let the point A 3 A U n B[ be denoted by E = I 3 (see Fig. 4.b). 

4. We set out also from the B\ ball configuration and we expand the horoball 
Pi until they comes into contact with their adjacent facets regarding P 24 
while keeping their 1 and 2-neighbouring horoballs tangent to them. More¬ 
over, we ’’blow up” the 3-neighbouring horoballs of Pi while their 1-neigh- 
bouring horoballs touch them. At this configuration e.g. the horoball Pi 
touches its neighbouring facet A 3 A 4 A 7 A 8 A n A 2 4 at the point T (see Fig. 3, 
and Fig. 4.b) and touch its 1-neighbouring horoballs e.g. P 3 ,P 5 ,P 7 ,P 9 
and its 2-neighbouring horoballs e.g. Pn. Furthermore, the ’’expanded” 
horoballs e.g. P 4 , P 6 , P 8 , Pi 0 touch the ’’shrunk” horoballs Pn and Pi 3 . 

This horoball arrangement is denoted by B 3 . 

5. Now we start from the configuration B 0 and we choose three arbitrary, mu¬ 
tually 3-neighbouring horoballs and expand them until they comes into con¬ 
tact with each other while keeping their 1 -neighbouring horoballs tangent 
to them. We note here that this horoball configuration can be realized in 
H 4 5 (see the subsection 4.2.4). At this configuration which is denoted by 
P 4 the horoballs are included on 2 classes related to P 24 , e.g. the horoballs 
Pi, P 10 , P 17 are in same type touching each other and their ’’smaller” 1- 
neighbouring horoballs that are also in same type. 
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4.1 Optimal horoball packings with horoballs in same type 

In this Section we consider the packings of horoballs where Vol(Bi D P 24 ) = 
Vol(Bj D P 24 ) for all i,j e {1,2,..., 24} thus the horoballs B, are in the same 
type regarding P 24 . 

It is clear that in this case the maximal density can be achieved if the neigh¬ 
bouring horoballs touch each other at the centres of the edges of P 24 and the 
density of this densest packing P 0 is equal to the maximal density of the horoball 
packings related to the Coxeter simplex tiling {3,4, 3,4}. For example in this case 
two horoballs P 4 and P 3 touch at the ’’midpoint” T\ of edge AiA 3 as projection 
of the polyhedron centre on it (see Fig. 3). These ball packings were investigated 
by the author in lf27l : 


V 0 := Vol(Bi n P 24 ) = -^\/2sinh ^arcosh^^jj « 0.00694, 


Vol(B 24 ) = 


7r 


864’ 


S(B 0 ) = 


Voi(Bi n P 24 ) 


(4.1) 


Vol{B. 


0.60793. 


24 J 


4.2 Optimal horoball packings with horoballs in different types 

The type of a horoball is allowed to expand until either the horoball comes into 
contact with other horoballs or with a adjacent facet of the honeycomb. These 
conditions are satisfactory to ensure that the balls form a non-overlapping horoball 
arrangement, as such the collection of all horoballs is a well defined packing in 

e 4 . 


4.2.1 Horoball packings B\ and their densities between the horoball ar¬ 
rangements Bo and B\ 

We set out from the Bo ball configuration (see above Section) and consider two 
1-neighbouring horoballs e.g. B\ and B 3 from it. Let J 0 = 1(0) — Tj be their 
point of tangency on side AiA 3 (see Fig. 3 and Fig. 4.a). Moreover, consider the 
point I(x) on the segment A 3 A 3 where the modified horoballs B t (x), (i = 1,3) 
are tangent to each other and x is the hyperbolic distance between 7(0) and I(x) 
(the value of x can also be negative if / (x) is on the segment TiA{). 

We blow up the horoballs P 4 (0) and Pn(0) (and also the horoballs P 2 , P 42 , 
B 14 , P 23 , P 24 and P 13 to achieve the B\ horoball configuration) until they come 
into contact with each other at the centre T 3 of octahedral facet AiA 3 A 3 A 7 A 9 Ah. 
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At this situation (see Fig. 3) the horoball centered at A\ is denoted, by -Bi(pi) 
where pi is the hyperbolic distance between I 0 and ii (see Fig. 4.a). 

The foot-point of the perpendicular from T 3 onto the staight line A 3 A 3 is J 0 = 
Ti which is the common point of the horoballs i?i(0) G £> 0 and f? 3 (0) G B 0 
centered at A ] and A 3 , respectively. The hyperbolic distance s 3 = T , T 3 between 
the points Ti[ti] and T 3 [t 3 ] can be computed by the formula (2.5) (see Fig. 4.a): 
The parallel distance of the angle 0 1 = T , T 3 A | Z is therefore we obtain by the 


T n =A, T 0 = A, 



Figure 4: Computations of hyperbolic distances p\ = I 0 I i and p 2 = I>T = I 3 T 3 

classical formula of J. Bolyai and by formula (2.5) the following equation (see 
Fig. 4.a). 

-= cosh si = \/2. (4.2) 

Sill (0i) 

We consider two horocycles "Ho and Hi through the points J 0 and I } with 
center A x in the plane /I i A 3 T 3 and the point Hi D A{T 3 is denoted by M. The 
horocyclic distances between points J 0 , M and I \, T 3 are denoted by h 0 and hi. 
By means of formula of J. Bolyai and of (4.2), we have 

^ = e pl = . . =>px = log(v/2) ~ 0.34657. (4.3) 

h 0 sm (0i) 

We extend the above modifications and denotations for all horoballs of packings 
between horoball arrangements B 0 and B\ i.e. the horoballs are denoted by B t {x) 
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(»e [o,Pi]). if x = 0 then we get the £> 0 horoball packing and if x — p\ then the 


£>i one. 

We obtain using the results of the former computations and of Lemma 4.2 the 
next 


Lemma 4.3 The density of packings B ( \ (see Fig. 5.a) between the main horoball 
arrangements Bq and B\ can be computed by the formula 


S(Bh(x)) 


TZ 0 Voi(Bfx) np 24 ) 
Vol(P 2A ) 


384 • V 0 (e 3x + 2 • e~ 3x ) 


x G [0,pi] 


and the maxima of function 8 (Bq(x)) (see Fig. 5.a) are realized at x — pi ~ 
0.34657 where the horoball packing density is S(Bl(pi)) ~ 0.71645. 




a. b. 

Figure 5: The graphs of functions S(Bq(x)) and S(Bf(x)) where x G [0, pf. 


Remark 4.4 We note, here that the above optimal density d{B' {] (f)\)) ~ 0.71645 
is equal to the density of known densest ball and horoball packings in H 4 . 

4.2.2 Horoball packings B\ and their densities between the horoball ar¬ 
rangements Bi and B 2 

We start our investigation from the B\ ball configuration. Here e.g. the horoballs 
B\ and B 3 touch each other at the point I\ (see Fig. 4.a) and Bi touch B n at 
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the point T 3 (see Fig. 3 and Fig. 4.b) etc. Furthermore, in B\ the common point 
of horosphere B{ with the line segment AfT is denoted by J 2 = /*(0) = D 
(see Fig. 4.b). We consider the point I*(x ) on the segment AfT where a modified 
horosphere B[(x) intersects the line segment AfT and x is the hyperbolic distance 
between P(0) and I*(x) (the value of x can also be negative if I* (x) is on the 
segment A 1 /*(0)). Corresponding to the above notions we introduce the notations 
# 13 (0) and B 13 (x). 

We blow up the horoballs B\ (0) and B \ 3 (0) while keeping their 1-neighbouring 
horoballs tangent to them until they comes into contact with their adjacent facets 
of P 2 4 e.g. upto the horoball B\ (x) touches the octahedral facet A 3 A 4 A 7 A 8 AnA 24 . 
At this arrangement relating to Fig. 4.b the horoball centered at /l, is denoted, by 
B 1 (p 2 ) where p 2 is the hyperbolic distance between J*(0) and T. 

The foot-point of the perpendicular from T onto the staight line A 3 A U is 
T 3 . The hyperbolic distance s 2 = TT 3 between the point T[t] and T 3 [ti] can 
be computed by the formula (2.5) (see Fig. 4.b): The parallel distance of the angle 
cp 2 = A\TT 3 Z. is s -2 therefore we obtain by the classical formula of J. Bolyai and 
by formula (2.5) the following equation (see Fig. 4.b): 

. ] , - = cosh s 2 = V2. (4.4) 

sm (<p 2 ) 

We consider two horocycles FL 2 and T~i 3 through the points I 2 and T with 
center A 3 in the plane AfTT 3 and the point 'H,> fl AfT 3 is denoted by E — I 3 . The 
horocyclic distances between points I 2 , T 3 and T, E are denoted by h 2 and h 3 . 
Similarly to (4.3) we obtain that p 2 = log(v / 2) ~ 0.34657. 

We extend the above modifications and denotations for all horoballs of pack¬ 
ings between horoball arrangements £>, and B 2 i.e. the horoballs are denoted by 
Bi(x) (i G [0, p 2 \). If x = 0 then we get the B\ horoball packing and if x = p 2 
then the B 2 one. 

We obtain using the results of the former computations and of Lemma 4.2 the 
next 


Lemma 4.5 The density of packings B\ (see Fig. 5.b) between the main horoball 
arrangements B\ and B 2 can be computed by the formula 


TZ 0 Vol(Bi(x) nP 24 ) 


5(Bl(x)) = 


Vol(P 24 ) 

48 • V 0 (2 e 3{pl+x) + 6 • e -^~ pl+x) + 16 • e ~^ pl+x) ) 


-7T 2 
3 ' 


x G [ 0 ,p 2 ] 
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and the maxima of function 5 (Bfx)) (see Fig. 5.b) are realized at x = 0 i.e. at 
the B{ ball packing (see Lemma 4.3). 

Remark 4.6 The density 5(Bfp 2 )) is equal to the maximal density of packing 
with horoballs in same types: 5(Bfp 2 )) = 5(Bf ~ 0.60793. 


4.2.3 Horoball packings B\ and their densities between the horoball ar¬ 
rangements Bi and B 3 

Similarly to the above subsection we set out from the B\ ball configuration and 
we will use the notations of subsection 4.2.2. Now, we expand the horoball -£>i(0) 
until they comes into contact with their adjacent facets regarding P 2 4 while keep¬ 
ing their 1 and 2-neighbouring horoballs tangent to them. Moreover, we ’’blow 
up” the 3-neighbouring horoballs of If (0) while their 1-neighbouring horoballs 
touch them. At this procedure this horoball is denoted by If (x). If we achieved 
the endpoint of this extension then e.g. the horoball Bfp 2 ) touches its neigh¬ 
bouring facet A 3 A 4 A 7 A 3 AnA 24 at the point T (see Fig. 3, and Fig. 4) and touch 
its 1-neighbouring horoballs e.g. B 3 , If,. If, B 9 and its 2-neighbouring horoballs 
e.g. B u . Furthermore, the ’’expanded” horoballs e.g. If, If,. If, If {) touch the 
’’shrunk” horoballs B u and B 13 . 

We extend the above modifications and notations for all horoballs of packings 
between horoball arrangements B\ and B 3 i.e. the horoballs are denoted by Bfx) 
(i e [0 ,P2\). If x = 0 then we get the B\ horoball packing and if x = p 2 then the 
B 3 one. Finally, we obtain the next 

Lemma 4.7 The density of packings Bf between the mean horoball arrangements 
B\ and B 3 can be computed by the formula 


m = = 

V ol(P 24 ) 

48 • Vo (e^ pl+ A + 7 • e -3(-pi+A _)_ 3 . e -3(pi+x) 3 . e -^{pi-x)\ 

-7T 2 

3 


, x e [ 0 , p 2 \ 


and the maxima of function 8 (Bfx)) are realized at x = 0 i.e. at the B\ ball 
packing (see Lemma 4.3). 


Remark 4.8 The function <5 (Bfx)) is the same with 5 (Bfx)) (x G [0, p 2 ]) (see 
Fig. 5.b). 
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4.2.4 Horoball packings £>d and their densities between the horoball ar¬ 
rangements Bq and B\ 

Here we consider the horoball configuration Bo and we choose three arbitrary, 
mutually 3-neighbouring horoballs e.g. B 4 , B w and £> 17 and let J 6 = /*(0) be the 
point of intersection of horosphere B{( 0) with the segment TA\ . Moreover, con¬ 
sider the point iAx) on the segment IfT where the expanded horosphere B{(x) 
intersects the segment IfT and x is the hyperbolic distance between /*(0) and 
I t (x) (see Fig. 6.a). We have seen in former subsections that the hyperbolic dis¬ 
tance between J 0 and T is 2pi = 2p 2 (see Fig. 5a and Fig. 5b). We consider a 
horocycles 'Hr, through the point T with center Ai in the plane A ] A ]i) T and the 
point "H .5 fl A 4 A 10 is denoted by K = I 5 . 

The foot-point of the perpendicular from T onto the staight line A 4 A W is 
called by Q whose coordinates are Q (l, ^=, 0, ^ ). 

We obtain by in the subsections 4.2.1 and 4.2.2 described method that the 
hyperbolic distance p 3 of the points 0 and K is p 3 = log y. 

The centre (’’midpoint”) of segment A 4 A 10 is denoted by 11 (see Fig. 6.a) (in 
our model this is Euclidean midpoint of segment AiA w , as well) whose distance 
p 4 to Q can be computed by the formula (2.5): p 4 = arccosh^-j^j . The point H 
lie on the line segment QK because 0.60199 ~ p 3 > p 4 ~ 0.45815. 

Finally, we obtain using the results of the former computations and of Lemma 
4.2 the next 


Lemma 4.9 The density of packings Bq (see Fig. 5.b) between the main horoball 
arrangements Bo and B 4 can be computed by the formula 


Y:t 0 yoi(B i (x)np 24 ) 


5(B*(x)) = 

48 • l/ 0 (3 • e 3x + 21 • e" 3 * 

-7T 2 

.3 


Vol(P 24 ) 

X e [0, 2pi + p 4 - p 3 ~ 0.54931] 


and the maxima of function S(Bq(x)) (see Fig. 6 .b) are realized at x = 0 where 
the horoball packing density is S(Bq(0 )) ~ 0.60793. 


Remark 4.10 The density S(Bq(0)) is equal to the maximal density of packings 
with horoballs in same types: 5 (Bq(p 2 )) = 6 (Bo) ~ 0.60793. 
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Figure 6: a. The computation of hyperbolic distance p 3 = IfT and p 4 = 
QH. b. The graph of function 8 (B^(x)) x G [0, 2pi + p 3 — p 4 ], 

4.3 Optimal horoball packings to hyperbolic 24-cell 

The main result of this paper is summarized in the following 

Theorem 4.11 The horoball arrangement B\ (see 4.2.1) provide the maximal 
horoball packing density related to the hyperbolic tiling V 24 with Schldfli sym¬ 
bol {3,4, 3,4} and its density is 8 opt (B) ~ 0.71645 if horoballs of different types 
are allowed at each asymptotic vertex of the tiling. 

Remark 4.12 The optimal horoball packing described and determined in this pa¬ 
per is a new horoball configuration which provide the known maximal density of 
realizable packings of the entire hyperbolic space H 4 . 


Proof 


It is well known that a packing is optimal, then it is locally stable i.e. each ball 
is fixed by the other ones so that no ball of packing can be moved alone without 
overlapping another ball of the given ball packing. 

The packings of horoballs can be easily classified by the type of ’’maximally 
large” horoball regarding the horoball packing to V-u- If we fix the ’’maximally 
large” horoball related to the above tiling then all possible horoball packing can 
be modified to achieve one of the above horoball configurations B] (x) ( 1 , j G 
{0,1, 2, 3,4}, i < j) without decrease of the packing density. 
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A horoball Bfx) is ’’maximally large” if Vol(Bi(x ) D P 24 ) (j £ 1... 24) is 
maximal. Here the maximal volume is denoted by Vol{B] nax ). 

1 . if ±voi(Br*) < V 0 then the maximal density can be computed by Sect. 
4.1 where the maximal density is 6 (Bo) ~ 0.60793. 

2. If V 0 < jg-Vol(B™ ax ) < V 0 -e 3pi then the optimal density can be computed 
by Sections 4.2.1, here the optimal density is 5(Bl(p 1 )) ~ 0.71645. 

3. If V 0 ■ e 3pl < P • Vol(B™' ax ) < Vq ■ e 6pi then the densities can be computed 
by Sections 4.2.2, 4.2.3 and 4.2.4 where the maximal density is 5(Bo) ~ 
0.60793. 

The volume of the ’’largest horoball” Vol(Bf ax ) < Vo • e 6pi therefore we 
proved the above Theorem. □ 

The above results also show, that the discussion of the densest horoball pack¬ 
ings and coverings in the n-dimcnsionalcn hyperbolic space with horoballs of 
different types has not been settled yet. Similarly to these, the problems of the 
densest hypersphere (or hyperball) packings and coverings are open, as well. 

Optimal sphere packings in other homogeneous Thurston geometries form 
also a class of open mathematical problems (see ll28l . 11291 , lf30l . OTft . If32l . If33ll . 
Ifl6l . [fTTl . G81). Detailed studies are the objective of ongoing research. 


References 

[i] 

Bezdek, K. Sphere Packings Revisited, European Journal of Com¬ 
binatorics (2006)27/6,864-883. 

[2] 

Bowen, L. - Radin, C. Optimally Dense Packings of Hyperbolic 
Space, Geometriae Dedicata, (2004) 104,37-59. 

[3] 

Bohm, J. - Hertel, E. Polyedergeometrie in n-dimensionalen 
Rdumen konstanter Kriimmung, Birkhauser, Basel (1981). 

[4] 

Boroczky, K. Packing of spheres in spaces of constant curvature, 
Acta Math. Acad. Sci. Hungar., (1978) 32,243-261. 

[5] 

Boroczky, K. - Florian, A. Uber die dichteste Kugelpackung im 
hyperbolischen Raum, Acta Math. Acad. Sci. Hungar., (1964) 15 , 
237-245. 


22 


J. Szirmai 


[6] Coxeter, H. S. M. Regular honeycombs in hyperbolic space. Pro¬ 
ceedings of the International Congress of Mathematicians, Amster¬ 
dam, , (1954) III, 155-169. 

[7] Fejes Toth, G. - Kuperberg, G. - Kuperberg, W. Highly Saturated 
Packings and Reduced Coverings, Monatshefte fur Mathematik, 
(1998) 125/2, 127-145. 

[8] Johnson, N. W. - Kellerhals, R. - Ratcliffe, J. G. and Tschantz, S. T. 
The size of a hyperbolic Coxeter simplex, Transformation Groups 
(1999) 4/4,329-353. 

[9] Kellerhals, R. Ball packings in spaces of constant curvature and 
the simplicial density function. Journal fiir reine und angewandte 
Mathematik , (1998) 494,189-203. 

[10] Kellerhals, R. Regular simplices and lower volume bounds for hy¬ 
perbolic n-manifolds, Ann. Globed Anal. Geom., (1995) 13,377- 
392. 

[11] Kolpakov, A. On the optimality of the ideal right-angled 24 cell, 
Algebr. Geom. Topol., (2012) 12/4,1941-1960. 

[12] Kozma, T. R. - Szirmai, J. Optimally dense packings for fully 
asymptotic Coxeter tilings by horoballs of different types, Monat¬ 
shefte fiir Mathematik, (2012) 168,27-47, DOI: 10.1007/s00605- 
012-0393-x. 

[13] Kozma, T. R. - Szirmai, J. New Lower Bound for the Optimal Ball 
Packing Density of Hyperbolic 4-space, Discrete Comput. Geom., 
(2015) 53, 182-198, DOI: 10.1007/s00454-014-9634-l. 

[14] Molnar, E. The projective interpretation of the eight 3-dimensional 
homogeneous geometries, Beitr. Algebra Geom., (1997) 38/2, 261— 
288. 

[15] Molnar, E. - Prok, I. - Szirmai, J. Classification of tile-transitive 
3-simplex tilings and their realizations in homogeneous spaces. 
A. Prekopa and E. Molnar, (eds.). Non-Euclidean Geometries, 
Janos Bolyai Memorial Volume, Mathematics and Its Applications, 
Springer (2006) Vol. 581, 321-363. 



Horoball packings related to hyperbolic 24 cell 


23 


[16] Molnar, E. - Szirmai, J. Volumes and geodesic ball packings to the 
regular prism tilings in SL 2 R space. Publ. Math. Debrecen, 84/1-2 
(2014), 189-203, DOI: 10.5486/PMD.2014.5832. 

[17] Molnar, E. - Szirmai, J. - Vesnin, A. Packings by translation balls in 
SL^R. J. Geometry 105/2 (2014), 287-306, DOI: 10.1007/s00022- 
013-0207-x. 

[18] Molnar, E. - Szirmai, J. - Vesnin, A. Geodesic and translation ball 
packings generated by prismatic tesselations of the universal cover 
of SL 2 (R). Submitted Manuscript 

[19] Marshall, T. H. Asymptotic Volume Formulae and Hyperbolic Ball 
Packing, Annates Academice Scientiarum Fennicce: Mathematica, 
(1999) 24,31-43. 

[20] Milnor, J. Geometry, Collected papers, Publish or Perish, (1994) 
Vol 1. 

[21] Radin, C. The symmetry of optimally dense packings, A. Prekopa 
and E. Molnar, (eds.). Non-Euclidean Geometries, Janos Bolyai 
Memorial Volume, Mathematics and Its Applications, Springer 
(2006) Vol. 581, 197-207. 

[22] Ratcliffe, J. G. and Tschantz, S. T. The volume spectrum of hyper¬ 
bolic 4-manifolds, Experiment. Math. (2000) 9/1, 101-125. 

[23] Rogers, C. A. Packing and covering, Cambridge University Press, 
(1964). 

[24] Slavich, L. Some hyperbolic 4-manifolds with low volume and 
number of cusps, Manuscript (2014) http://arxiv.org/abs/1402.2580 

[25] Szirmai, J. Horoball packings for the Lambert-cube tilings in the 
hyperbolic 3-space, Beitr. Algebra Geom., (2005) 46/1, 43-60. 

[26] Szirmai, J. The optimal ball and horoball packings of the Coxeter 
tilings in the hyperbolic 3-space, Beitr. Algebra Geom., (2005) 46/2, 
545-558. 


24 


J. Szirmai 


[27] Szirmai, J. The optimal ball and horoball packings to the Cox- 
eter honeycombs in the hyperbolic ("/-space. Beitr. Algebra Geom., 
(2007) 48/1, 35-47. 

[28] Szirmai, J. The densest geodesic ball packing by a type of Nil lat¬ 
tices, Beitr. Algebra Geom., (2007) 48/2, 383-397. 

[29] Szirmai, J. The densest translation ball packing by fundamental lat¬ 
tices in Sol space, Beitr. Algebra Geom., (2010) 51/2, 353-373. 

[30] Szirmai, J. Geodesic ball packing in S 2 x R space for generalized 
Coxeter space groups, Beitr. Algebra Geom., (2011) 52, 413-430. 

[31] Szirmai, J. Geodesic ball packing in H 2 x R space for generalized 
Coxeter space groups. Math. Commun., 17/1 (2012), 151-170. 

[32] Szirmai, J. Lattice-like translation ball packings in Nil 
space. Publ. Math. Debrecen, 80/3-4 (2012), 427-440, DOI: 
10.5486/PMD.2012.5117. 

[33] Szirmai, J. A candidate to the densest packing with equal balls in 
the Thurston geometries. Beitr. Algebra Geom., 55/2 (2014), 441- 
452, DOI 10.1007/sl3366-013-0158-2. 

[34] Szirmai, J. Horoball packings and their densities by generalized 
simplicial density function in the hyperbolic space, Acta Math. 
Hungar., (2012) 136/1-2, 39-55, DOI: 10.1007/sl0474-012-0205- 
8 . 

[35] Szirmai, J. Horoball packings to the totally asymptotic regular sim¬ 
plex in the hyperbolic n-space, Aequat. Math., 85 (2013), 471-482, 
DOI: 10.1007/s00010-012-0158-6. 



